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Abstract 

We show that the double asymptotic scaling of the HERA structure function data is consistent 
with pre-HERA data at larger x, soft pomeron behaviour at small x and a sensible starting scale 
Qo. We can thus actually calculate at small x and large Q 2 by evolving up perturbatively at 
two loops, without any fitting. 
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Perturbative QCD makes a definite prediction |l| for the shape of the rise of the 
singlet part of the structure function F 2 at large Q 2 and small x, provided only that the 
starting distribution is not too singular. This prediction is most striking when expressed 
as a double asymptotic scaling in the two scaling variables || 

where t = ln(Q 2 /A 2 ). When F 2 is rescaled by a multiplicative factor 

R F = Na l ' 2 pe- 2lC7+5a / p , (2) 

it should scale in both a and p as they become large, since the growth predicted in [|1| 
has been scaled out. The parameter 7 which controls the rate of growth of F 2 is simply 
related to the leading coefficient of the beta function (in fact 7 = ^/12//3q), while 5 is an 
anomalous dimension. Recent data from HERA || are in remarkably good agreement with 
this double scaling prediction, to the extent that it may be used to directly determine /3q 
by measuring the slope of the riseQ. It thus now becomes imperative to compute higher 
order corrections: double scaling violations. 

The original calculation in [fy expanded the one loop anomalous dimensions around 
their leading singularity, which should dominate the evolution at large Q 2 and small x 
if there is no corresponding singularity in the input distribution. In [[2] we showed that 
this approximation turns the Altarelli-Parisi equation for the gluon distribution xg(x; t) 
into a two-dimensional wave equation with light-cone variables ln(xo/x) and \n(t/to). The 
(unstable) propagation of soft boundary conditions (in particular g(x;to) ~ x _1 as x — > 0, 
the intercept of the soft pomeron being close to unity) into the interior of the light-cone 
then produces the generic rise in F 2 . Asymptotically the details of the boundary conditions 
are washed out, and only their overall normalization is left. The simplest approximation is 
thus sufficient to capture the essential physics at small x: F 2 rises because of the instabilty 
of gluons to gluon emission via the triple-gluon vertex. 

Of course the success of this simple picture still leaves some important questions 
unanswered. In particular, is the double asymptotic scaling behaviour consistent with the 
structure functions measured (pre-HERA) at larger values of xl Furthermore, how good 
is the leading singularity approximation to the splitting function, and how large are the 
two loop corrections? In other words, it would be useful to compute the normalization and 
sub-asymptotic corrections to double scaling, not only to explain why double scaling works 
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so well, but also to refine the comparison with the increasingly precise experimental data. 
Post-asymptotic corrections due to singularities in yet higher orders of perturbation theory, 
or to higher twist corrections from parton recombination effects, are discussed briefly in 
H; they will not be considered here. 

To compute the normalization and sub-asymptotic corrections numerically, we adopt 
the following four step algorithm: 

(a) Take a set of parton distributions A which has been fitted to pre-HERA data 
with Q 2 £ 4GeV 2 , x £ 10" 2 . For illustration we will use here the MRS SO', DO' and D-' 
distributions @. 

(b) Evolve A backwards to a starting scale Qq. This scale should be chosen sufficiently 
low that it is reasonable to match to soft non-perturbative behaviour there, but not so low 
that perturbation theory has become meaningless. Here we choose Qq = lGeV, as in 0. 

(c) Remove the (unphysical) small-x tails of the distributions (that is those parts with 
x ^ 10 ~ 2 , and replace them the conventional expectation from Regge theory. In particular 
the glue and sea distributions should be given soft tails, xg(x,to) ~ a; -0 08 , to match the 
intercept of the (nonperturbative) soft pomeron. This gives new distributions, which we 
will call (A)0. 

(d) The new distributions (A)0 are now evolved back up to high Q 2 , where they can 
be compared to experimental data (at low Q 2 such comparison would be difficult because 
of contamination by higher twist corrections). At large-x they will (by construction) be 
indistinguishable from the original distributions A. At small- x they should now exhibit 
double scaling, but with a definite normalization and sub-asymptotic corrections. They 
can thus be compared directly to the HERA data. 

The results of this procedure^ are shown in Fig. 1 (the distribution (D0')0; (S0')0 is 
very similar) and Fig. 2 (the distribution (D-')0), presented as a and p scaling plots; they 
should be compared with the scaling plots in .@ The original MRS distributions DO' 
and D-' are also shown (dotted; see also ||). The figures largely speak for themselves. 
The normalization of the HERA data || is correctly reproduced, with no free parameters 
(unless Qq were to be considered as such): this is not a fit to the HERA data! The 
subleading and two loop corrections do not spoil double scaling: although the slope of the 

1 The parameters of the three new sets of distributions functions (S0')0, (D0')0 and (D-')0 
may be obtained by email from rball@suryall . cern. ch. 

2 Note that in the p plots most of the data points now have a 3> 1.1. 
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x 2 


Best fit Nc 

x 2 


trmalization 
Rel. Norm. 


% mom. 
in glue 


DAS§ 




107 




33% 


(SO')O 
(D0')0 
(D-')O 


115 [309] 
126 [263] 
114 [323] 


105 [230] 

106 [213] 
112 [185] 


1.04 [1.13] 
1.06 [1.10] 
0.98 [0.87] 


40% [44%] 
41% [44%] 
38% [43%] 



Table. 



a-plots is now a little lower (2.09 ± 0.35 for (D0')0, 1.92 ± 0.35 for (D-')0) than at leading 
order (where it was 2.4; the data have a slope of 2.37 ± 0.16), there is now a slight rise in 
the p plot at large p. 

To show quantitatively how well the new distributions account for the data, we also 
give a table of values of x 2 (f° r all 115 data points with x < 0.1). In the third column we 
show how the \ 2 f ans if the normalization of the distributions is left free. For comparison 
we give in square brackets similar statistics for the original MRS distributions ||, none of 
which fit the HERA data because they do not exhibit double scaling (SO' and DO' because 
Qo was too large, D-' because it incorporates the power-like singular growth inspired by 
the Lipatov 'hard pomeron').! The more recent H and A distributions @, which achieve a 
fit to the HERA data by introducing two new free parameters (basically parameterizing a 
small admixture of a power rise at small- a;, thus roughly interpolating between DO'and D-') 
have for comparison a x 2 °f 100- From the calculations presented here it should be clear 
that similar (perhaps better) results could have been obtained with no new parameters 
simply by dropping Qq from 4GeV 2 to lGeV 2 , and there taking a soft initial distribution 
of the same form as DO': indeed if such an approach had been originally taken in ||, an 
extremely accurate prediction for Fl± at small x would have resulted. 

If using perturbation theory at scales as low as lGeV makes one uncomfortable, one 
could instead incorporate the expected asymptotic behaviour |j] by hand into a starting 
distribution at 4GeV 2 ; it really does not matter how the double scaling is produced, 
provided it obtains throughout the HERA kinematic region. After all, the most striking 
feature of the small-x data is the precocious onset of double asymptotic scaling 0. But 
if one evolves at two loops from a soft distribution at lGeV, one can also generate both 
the correct normalization and subasymptotic corrections. 

3 The GRV distributions Q also fail; they implicitly incorparate double scaling, but overshoot 
the data because their starting scale is much too low. 
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Figure Captions 



Fig. 1. Double scaling plots of RfF^ vs. a) a and b) p. The data are taken from ref.0, 
and the curves are those of the new parton distributions (D0')0. The old MRS 
prediction DO' is shown dotted. 

Fig. 2. As fig. 1, but with the curves now corresponding to (D-')O and D-'(dotted). 



6 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409374vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409374vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409374vl 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409374vl 
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